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In this Letter, we consider a neutral system of mobile positive and negative charges confined on 
the surface of curved films. This may be an appropriate model for: i) a highly charged membrane 
whose counterions are confined to a sheath near its surface; ii) a membrane composed of an equimolar 
mixture of anionic and cationic surfactants in aqueous solution. We find that the charge fluctuations 
on ■ contribute a non-analytic term to the bending rigidity that varies logarithmically with the radius 

of curvature. This may lead to spontaneous vesicle formation, which is indeed observed in similar 
£^ ' systems. 
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fSJ ' Electrostatics of charged objects such as polyelectrolytes and membranes in aqueous solution plays an important 
role in many biological systems jj] . The fundamental description of these systems has been the mean-field approaches 
- the Poisson-Boltzmann (PB) or Debye-Hiickel (DH) theory (for a review, see Q). However, for a highly charged 
surface, the Manning theory of counterion condensation jj| provides an analytically tractable approximation to the PB 
theory. Indeed, it has been demonstrated rigorously from the solutions to the PB equation Q that the electrostatic 
potential far away from the charged surface is independent of the charge density above a certain critical value, implying 
that the counterions are confined to a thin layer close to the charged surface. However, like the PB theory, it fails 
Ch to capture the correlation effects of the counterions since it expressly assumes that the "condensed" counterions are 
j uniformly distributed. On physical grounds, we should expect that at low enough temperatures the fluctuations of 
<**] ■ these condensed counterions about a uniform density would give rise to new phenomena. Indeed, recent simulations 
||[| show that the effective force between two like-charged rods and planar surfaces actually becomes attractive at 

1 short distances. These suprising results shed new light on the understanding of the electrostatic adhesion between 
cells [Q and the puzzling problem of DNA condensation 0. In this Letter, we examine the effect of fluctuations of 

s ! i these condensed counterions on the bending rigidity of a charged membrane. 

The elastic properties of a fluid membrane are characterized by three macroscopic parameters - a bending elastic 
modulus k, a Gaussian modulus Kg, and a spontaneous curvature H . The deformation free energy per unit area, 
' expressed in terms of the mean curvature H and Gaussian curvature K may be given by the Helfrich free energy 

O ' 

f=^(H H Q f + K G K. (1) 

Within an additive constant, the free energy of a sphere with radius R is given by f s — (2 k + kc)/R 2 — 2kH /R and 
of a cylinder with radius R by f c — n/2R 2 — kHq/R. Therefore, the parameters kHo and k + k g may be determined 
I ■ from f s and f c . 

The problem of the electrostatic contribution to the bending constants of layered membranes within the PB mean 
field approach has been studied [ pT| . The electrostatic renormalization of the bending rigidity turns out to be 
positive; hence electrostatics augments the rigidity of charged membranes. Here we go beyond these PB approaches 
by assuming that the surface charge density n is sufficiently high that the condensed counterions are confined to a 
layer of thickness A << L, where A is the Gouy-Chapman length, which scales inversely with uq and L is the linear 
size of the charged membrane. By considering the in-plane fluctuations of the condensed counterions and charges 
on the membranes, we model the system effectively as a 2-D coulomb gas interacting with a r -1 potential. This 
model has yet another experimental realization - a neutral membrane composed of a dilute mixture of anionic (-) 
and cationic (+) surfactants. 

The electrostatic free energy of the system is the sum of the entropy of the charges and the electrostatic interaction 
energy among them: 
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where At is the de Broglie wavelength of the charges, Ib = € £ bT ~ 7 A is the Bjerrum length for an aqueous solution 
of dielectric constant e = 80 (H2O), /3 _1 = ksT, ks is the Boltzmann constant, T is the temperature, and n.;(x) is 
the coarse-grained density of the charges of species i. The domain of the integral in Eq. (^|) spans the entire charged 
membrane. In order to calculate the change in the free energy due to fluctuations, we assume that ?ii(x) = ?iq + 5ni(~x) 
and expand the electrostatic free energy to second order in Srii p^ ]: 

Ib , ^ 2 (x 



(3AF e 



d 2 xd 2 x' 



x') 



2 n 



(3) 



where 5a = 5n + — <5n_. The first term in the bracket is the Coulomb interaction of the charges. The second term 
comes from the second variation of the ideal gas entropy of the charges. The change in the free energy is obtained by 
summing all fluctuations weighted by the Boltzmann factor: 



(3G e 



In 



J d5a(x) exp —/3AF 



(4) 



It should be mentioned that Eq. (Q) contains a divergent self-energy term which has to be substracted out. This 
means that we have to discard the first two terms in the expansion for Is 0, as can be seen easily by considering 
the zero temperature limit. As T — > 0, the free energy is reduced to the electrostatic energy which is first order in 
Ib- Since the self-energy is just a constant independent of temperature, it must be linear in Ib- In the following, we 
employ this "substraction scheme" together with Eq. (Q) to calculate the free energy where charges are confined to 
the surfaces of three geometries: i) a plane, ii) a sphere, and iii) a cylinder. 

For the case of charges confined to a plane AF in Eq. (|3|) can be diagonalized by Fourier transform and is quadratic 
in <5(7. Performing the Gaussian integrals in Eq. (^) and substracting out the self-energy term, we obtain the free 
energy per unit area due to fluctuations |p^ , |l3| 
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where (3 — l/fc^T and 1/Xo = Sirnols, which scales like the Gouy-Chapman length, is a length scale analogous to 
the Debye screening length in 3-D. This result can also be obtained by solving the Debye-Hiickel equation in 2D [fLjj . 
Note that Eq. (J5J) is ultravioletly divergent because of the infinite energy associated with the collapse of opposite 
charges. Thus a microscopic cut off is necessary. In ref. [ {l5| , the author shows by partial summation of the Mayer 
series that the resulting free energy is convergent and indeed equivalent to a microscopic cut-off. From Eq. (^), one 
can deduce that the screening of the charges is weak and the potential to the lowest order has the distance dependence 
of a dipolar field, where A^ plays the role of the dipole moment. 

For the case of charges confined on a sphere of radius R, after following a similar procedure to that described above 
we obtain 
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It is easy to show that by setting k — l/R and taking the limit R — > 00, we recover the planar result. Equivalcntly 
we may write Eq. (||) as 
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The difference g sp — g p i , can be evaluated as an asymptotic expansion in 1 / R using the Euler-MacLaurin summation 
formula @] with f(l) = (21 + 1) ln(2Z + 1 + R/X D ). The result is: 
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In deriving the result above, we have regularized the integral in Eq. (j?]) and the sum in Eq. (^|) by an ultraviolet 
cut-off A. However, the leading term in Eq. (||) is cut-off independent and those higher order cut-off dependent 
terms tend to zero as A — > 00. 

For the case of a cylinder, we obtain the free energy: 
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where I m and K m are modified Bessel functions of order m. The evaluation of the integrals here is relatively difficult. 
However, we argue that g cy i — g p i has the following asymptotic expansion: 

P(g cy i- g P i) = -4^2 HR/ad ) + .-, (10) 

for R — > 00. First, we note that the only relevent contributions to the g-integral in Eq. (^|) are sharply peaked at 
q « with width A q ss m/R. Hence, the Bessel functions can be approximated by I m (qR)K m (qR) ~ l/2m, yielding 
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Equation ( fTo|) can now be obtained by using the Euler-MacLaurin summation formula with f(m) = m log(2m+ R/ Xd) . 
The modificaitons to the bending constants can be obtained from Eqs. (||) and ( |l0| ) to yield 

Ak = -^ In(i2/A D ), 

Ak g = --i- ln(i?/A D ). (12) 
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We thus find that the contribution to the membrane elastic contants due to charge fluctuations is non-analytic. This 
kind of non-analyticity in the bending constants exists in the literature in other situations, for example in a system 
consisting of a membrane and rod-like cosurfactants |0. In the present case, this non-analyticity can be considered 
a signiture of 2-D charged systems. The DH theory in 3-D yields an expression for the change in the free energy per 
unit volume [fl8[ Af^D ~ — A^ 3 + .... In constrast, Eq. (^J) has a similar expansion for the free energy per unit area 
but contains a logarithmic term A/2D ~ — A^, 2 lnfA^/a] + .... Therefore, it is not unexpected to find logarithmic 

corrections to the bending constants. 

Secondly, we remark that both Ak and Akq are negative, in contrast to the mean-field PB contributions, where the 
renormalization of the bending moduli are always positive and the Gaussian moduli may be negative in some cases. 
In a system in which R/Xd » 1, An is large compared to the mean-field contribution and the membrane becomes 
more flexible. Therefore, charge fluctuations induce bending of a charged membrane. This conclusion can be seen 
physically in the light of the recent works Jl2],[l3| on attractive interactions between two planar charge-fluctuating 
membranes. It is found that for large distance h separating the two membranes, the attractive force per unit area 
scales as h~ 3 . Now, a sphere or a cylinder can be approximated as two flat planar surfaces in the limit R — > 00 and 
their interaction free energies per unit area therefore should scale like / ~ — R~ 2 . Hence negative contribution to the 
bending modulus. 

The negative contribution of Akq from charge fluctuations can also be understood physically as follows. Since the 
r _1 potential is rotationally invariant, one might expect that the charged surface has rotational symmetry when the 
the electrostatic free energy is a minimum. Now, recall that the Gaussian curvature, expressed in terms of the principal 
radii R\ and R2, is K = R ^ Ro ■ Therefore, as can be easily seen, the free energy in Eq. ([!]) is indeed lowest (since 
Aug < 0) if there is rotational symmetry about an axis normal to the charged surface or R\ = R2, in accordance 
with the rotational symmetry of the r _1 potential. Furthermore, the fact that Akq < has interesting experimental 
consequences since strongly negative values of kg favor the formation of many disconnected pieces with no rims, 
like spherical vesicles. Therefore, when the surface charge density is made sufficiently large, the membrane might 
spontaneously form vesicles, due to fluctuations of condensed counterions. Experiments [pL9|| on charged surfactant 
systems supports this conclusion. 

The result presented in this Letter is particularly relevant to recent experiments |20| where the authors find the 
formation of vesicles by mixing anionic and cationic surfactants. Two aspects of their experiment can be qualitatively 
accounted for by the present model. They find, in equilibrium, large vesicles with R ~ 1000 A and substantial size 
polydispersity. Indeed, the vesicle free energy per unit area given by 

f ves = Kb /R 2 - — ^-j ln( R/X D ), (13) 

where Kb is the bare value of the bending rigidity, has an equlibrium value R* ~ A^> e~xp(Kb/kBT), which can be 
large even for a moderate value of Kb of the order of 3-5 fcsT. Furthermore, the second derivative of the free energy 
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f"(R*) ~ e - K b/k B T j g ex p 0nen tially small. Hence the variance or fluctuations in R, < (AR) 2 >~ 1/ f"(R*) is large, 
implying size polydispersity. 

In conclusion, by studying fluctuations of charges on curved films, we have deduced non-analytic contributions 
to the bending energy of a membrane. Our calculation is applicable to condensed counterions on a highly charged 
membrane, and mixing of surfactants of opposite charges. For the latter case, we find qualitative agreements with 
experiments. 
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